The electro-disintegration of few body systems revisited 
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Recent studies of the electro-disintegration of the few body systems at JLab have revived the 
field. Not only recoil momentum distributions have been determined in a single shot. But also they 
confirm that the diagrammatic approach, which I developed 25 years ago, is relevant to analyze 
them, provided that the Nucleon-Nucleon scattering amplitude, determined in the same energy 
range, is used. They provide us with a solid starting point to address the issue of the propagation 
of exotic components of hadrons in nuclear matter 
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The primary goal of the study of the (e,e'p) reaction 
on nuclei was, and still is, the determination of the high 
momentum components of the nuclear wave function. In 
the past, the spectral functions measured at Saclay or 
Amsterdam suffered from large corrections (about a fac- 
tor two or more) due to Final State Interactions (FSI) 
and Meson Exchange Currents (MEC). A surveyof the 
state of the art at that time can be found in ref. [ij . The 
corresponding experiments were performed at low values 
(~ 0.4 GeV^) of the virtuality of the exchanged pho- 
ton. 

When it was decided to build CEBAF, a common be- 
lief was that increasing was the way to suppress FSI 
and MEC contributions. This is partly true, since both 
the FSI and MEC amplitudes involve a loop integral, 
which connects the nuclear bound and scattering states 
and which is expected to decrease when increases as 
form factors do. But this is partly wrong, since the sin- 
gular part of the FSI integral does not depend on , 
besides the trivial momentum dependency of the elemen- 
tary operators. It comes from unitarity, and corresponds 
to the propagation of an on-shell nucleon. It involves 
on-shell elementary matrix elements and it is maximum 
when the kinematics allows for rescattering on a nucleon 
at rest In the (e,e'p) channel, this happens in quasi- 
free kinematics, when X — jlrav ~ \ [y being the 
energy of the virtual photon, and m the nucleon mass). 

In turn, this kinematics provides us with a way to iso- 
late NN scattering (or more generally scattering between 
hadrons) and opens up an original use of the (e,e'p) reac- 
tions _3;, 4J : the study of exotic components of the hadron 
wave function via color transparency or color screening, 
for instance. 

Fig. n exhibits these features. It shows the angular 
distribution, against the neutron angle Qji with the vir- 
tual photon, of the ratio between the full cross section 
of the D(e,e'p)n reaction and the quasi-free contribu- 
tion, when the momentum Pr of the recoiling neutron 
is kept constant. FSI (dashed curves) are maximum 
near Qb. = 70° where X = 1 and on-shell rescattering 
is maximized. At low values of the recoil momentum 
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FIG. 1: The ratio between the full cross section and the con- 
tribution of the quasi-free scattering. 



{Pji — 200 MeV/c), on-shell nucleon rescattering reduces 
the quasi-free contribution, as expected from unitarity (a 
part of the strength of the quasi-elastic channel is trans- 
ferred to inelastic ones). At high values of the recoil mo- 
mentum (Pfj = 500 MeV/c) the quasi- free contribution 
strongly decreases as the nucleon momentum distribu- 
tion: on-shell rescattering takes over and dominates. 

Similarly, the A, which is produced on one nucleon and 
exchanges a meson with the second nucleon in the MEC 
amplitude, can also propagate on-shell. The correspond- 
ing singularity appears at larger recoil angles and shifts 
the NN rescattering peak (full curves) . In fact other bary- 
onic resonances can be excited and propagate, widening 
the peak further toward larger angles. But the A is the 
most prominent part of the nucleon response function, 
and the effects of the higher mass resonances are expected 
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to be smaller, except maybe at higher recoil momenta. 

Experimenls recently performed at JLab confirm 
this behavior, which was already predicted and mea- 
sured ^8] in the ttN rescattering sector at lower energy. 

To be more specific, the method ,2] is based on the 
expansion of the amplitude in terms of few relevant dia- 
grams, which are computed in the momentum space, in 
the Lab. frame. The kinematics as well as the prop- 
agators are relativistic and no angular approximation is 
made in the evaluation of the loop integrals. The elemen- 
tary operators which appear at each vertex have been 
calibrated against the corresponding channels. Its ap- 
plication to the D(e,e'p)n channel has been discussed in 
refs [1,^3 ^^'^ to ttis '^He(e,e'p) channels in refs. [lllll^ . 
A comprehensive summary is given in ref. for the 
'*He(e,e'p)T channel. 

For the sake of the discussion, I reproduce the Plane 
Wave (PW) and FSI amplitudes for the D(e,e'p)n channel 

Tpw = ^{mi\Jp{q^)\mp){^mp^m2\lMj)Uo (^2) — 1= 
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= 2.35 GeV 
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ref. |3| for a full discussion): This happens in the quasi- 
elastic kinematics, X = 1. The width of the on-shcU peak 
in Fig. 121 reflects the Fermi distribution of the target nu- 
cleon, while the ofF-shell (principal) part of the integral 
vanishes at X = 1. 
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where £'„ 



+C/2(| n \){2mi\m, \ lMj)r2""(i)}(2) 

\j m? + (fc — n)2 and p^ ~ Me, + v — 

\J im? + fp . The momenta and magnetic quantum num- 
bers of the outgoing proton and neutron are respectively 
Pi, p2, mi and m2, while the magnetic quantum num- 
ber of the target deuteron is Mj. The S and D parts 
of the deuteron wave function are respectively Uq and 
1/2- The relativistic expressions of the proton Jp{q^) and 
neutron Jn{q^) currents are used in both the PW and 
FSI amplitudes, contrary to where their expansion 
up to and including terms of order 1 /m^ was used: the 
difference does not exceed a few per cent, except at very 
forward or backward recoil angles. The FSI integral runs 
over the momentum n of the spectator nucleon. Since 
the energy is larger than the sum of the masses of the 
two nucleons, the knocked out nucleon (p, Ap) can prop- 
agate on-shell. Due to the dominance of the S'-wave part 
of the wave function, the corresponding singular part of 
the integral is maximum when the scattering of the elec- 
tron on a nucleon at rest is kinematically possible (see 



FIG. 2: The angular distribution of the on-shell (dot-dashed) 
and off-shell (dashed) parts of the ratio of the FSI to the PW 
cross sections. 



The physical picture is the following. The electron 
scatters on a proton at rest which propagates on-shell 
and rescatters on the neutron which is also at rest. In 
the Lab. frame, the soft neutron recoils at 90° with re- 
spect to the fast proton which is emitted in the forward 
direction. Two body kinematics imposes that the angle 
of the rescattering peak (dip) moves with the recoil neu- 
tron momentum: around 70° when pn = 500 MeV/c, 
80° when pn — 200 MeV/c. The same occurs, in a differ- 
ent part of the phase space, when the electron interacts 
with the neutron. In the classical Glauber approxima- 
tion, the nucleon propagator in Eq. El is linearized and 
recoil effects are neg lected: Therefore the rescattering 
peak stays at 90° '14, This drawback has been cured 
in the Generalized Eikonal Approximation (GEA) 
which takes into account higher order recoil terms in the 
nucleon propagator, and neglects only terms of the order 
p\/m?. It comes as no surprise that GEA predicts the 
FSI peak at the same place as in my diagrammatic ap- 
proach which takes into account the full kinematics from 
the beginning t2„ X3]. While it is valid at forward angles, 
the classical Glauber treatment is simply not correct for 
analyzing the (e,e'p) reactions at large angles and large 
recoil momentum. 

Since it involves on shell matrix elements and relies on 
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the low momentum components of the wave function, the 
FSI amplitude is founded on solid ground near X ~ 1, 
provided the correct parameterization of the NN ampli- 
tude is used. 

In the pre-CEBAF era, the relative kinetic energy of 
the two outgoing nucleons {Tl = Q^/2m ~ 200 MeV) 
was low enough to rely on the partial wave expansion 
of the nucleon-nucleon scattering amplitude (T^r a?), see 
for instance 0, , of which both the on-shell and half 
off-shell parts were solutions of the Lippman-Schwinger 
equation with the same potential (Paris) as for the bound 
state S, P and D waves were retained and the FSI 
loop integral was done analytically according to ref . ^20,] , 
fully taking into account Fermi motion effects (unfactor- 
ized calculation). When this is done, and the momenta 
expressed in the rest frame of the neutron-proton system, 
Eq. El coincides with Eg. C.8 of 
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FIG. 3: The variation of ctjvjv and /3 with energy. 



At higher energies (let's say when the relative kinetic 
energy of the outgoing fragments exceeds 500 MeV or 
so) , too many partial waves enter into the game and their 
growing inelasticities prevent us to compute the scatter- 
ing amplitude from a potential. It is better to use a 
global parameterization of the NN scattering amplitude. 
On general grounds it can be expanded as follows 



T-i 



NN 



a + 17(01 -t- (72) • k± + spin — spin terms (3) 



where k±_ is the unit vector perpendicular to the scatter- 
ing plane. 

Above 500 MeV, the central part a dominates. It is 



almost entirely absorptive, and takes the simple form 

a = — (e -I- i) aNN exp[-ij (4) 

In the forward direction its imaginary part is related to 
the total cross section a^N, while the slope parameter 
/3 is related to the angular distribution of NN scattering. 
Both can be determined from the experiments performed 
at Los Alamos, Saturne and COSY. Fig. O shows the 
values which I use. Below 500 MeV, I have extrapolated 
them in such a way the absorptive part of the amplitude 
vanishes at the pion production threshold. The ratio e, 
between the real and imaginary part of the amplitude, 
is small: I keep it constant (e = —0.2) above 1 GeV, 
and smoothly extrapolate it down to zero at the pion 
threshold. 

Such a parameterization is very convenient to compute 
the rescattering amplitude. It adds its absorptive part, 
which dominates at high energy, to its expansion in terms 
of the real part of phase shifts (of which I use the exper- 
imental values, above = 500 MeV), which dominates 
at low energy. However, at high energies, it leads only to 
an accurate prediction of its singular part (on-shell scat- 
tering). Contrary to low energy, there is unfortunately no 
way to constrain the half-off shell behavior of the absorp- 
tive part of the NN scattering amplitude, and one can get 
only an estimate of the principal part of the rescattering 
amplitude. It turns out that it vanishes at X=l (Fig. 
and it does not dominate at high energy. So, the method 
is founded on solid grounds in the quasi-elastic kinemat- 
ics (X^l). Away, it tells us in which kinematics FSI are 
minimized. 

The A formation and MFC amplitudes are computed 
according to eqs CI and C2 in ref 13J . I updated them by 
implementing the full relativistic expression of the ttNN 
vertex and using the latest dependency of the N ^ A 
electromagnetic form factor 2lj . Since it falls down more 
rapidly than the Nucleon form factor, the A formation 
amplitude is suppressed at high Q^. Also the unitary 
singularity associated with the A propagation is weaker 
than in the FSI amplitude since the A pole is distant 
from the energy axis by its half width. Again the A 
propagates almost on shell in the kinematics of Fig. ^ 
and it is worth to emphasize that the parameters are 
those which reproduce the NN A^A cross section in 
the few GeV range (see e.g. 22]). 

Fig. 0] shows the full angular distribution of the 
D(e,e'p)n reaction for = 5 GeV^, at the top of the 
unitary peak in Fig^ X ~ 1. The A formation term 
contributes little up to p„ ^ 800 MeV/c, but dominates 
above. At the extreme backward proton emission angles 
(large momentum of the neutron but vanishing momen- 
tum of the proton) the interaction of the electron with the 
neutron takes over and is modified, as the forward pro- 
ton peak, by FSI, MEC and A formation term. These 
findings are reproduced by the preliminary analysis of 
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FIG. 4: The momentum distribution in the D(e,e'p)n reaction 
at X=l and = 5 GeV^. Dashed hue: PW. Dash-dotted 
hue: with FSI. FuU hue: MEG and A included. 



the D(e,e'p)n reaction |3 recently recorded in the full 
phase space with CLAS at JLab. We must await its final 
analysis for a detailed comparison. 

So far, only the analysis of the electro-disintegration 
of "^He and ^He have been completed at JLab. Fig. |S1 
shows how well the diagrammatic method reproduces 
the cross section of the '^He(e,e'p)d reaction recently 
measured |3 with two magnetic spectrometers, at 
= 1.55 GeV^, in the quasi-free kinematics (X=l). 
The wave function is the solution 25l (2 Q1 of the Fad- 
deev equations for the Paris potential [ij. The nucleon 
single scattering (FSI) and two body MEG amplitudes 
are implemented as described in [13|. Both pp as well 
as T = and T — 1 np active pairs are considered. At 
such a high virtuality, the relative kinetic energy between 
the outgoing proton and deuteron is Tl — 830 MeV, 
where the NN cross section reaches its maximum and be- 
comes flat around unn — 45 mb. Again, FSI reduces 
the quasi-free contribution below 300 McV/c and over- 
whelms it by more than a factor five around 500 MeV/c. 
Above 1 GeV/c, MEG and A production enhance the 
cross section, but are unable to reproduce the last three 
experimental points around 1 GeV/c. Here, one en- 
ters into the kinematical regime where the deuteron is 
fast and emitted in the forward direction while the pro- 
ton is slow and becomes a spectator: this is responsi- 
ble for the small deuteron knockout peak at the extreme 
right of the figure. In order to accommodate the experi- 
ment around 1 GeV/c and above, one needs a mechanism 
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FIG. 5: The momentum distribution in the ^He(e,e'p)d reac- 
tion at X=l and = 1.55 GeV^. Dashed line: PW. Dotted 
line: with FSI. Dash-dotted line: 2 body MEG and A in- 
cluded. Full line: 3 body mechanisms included. 



which shares the photon momentum between the three 
nucleons. Three body meson rescatterings, computed as 
in ref. iii the right direction but fall short. It is 

very likely that nucleon double scattering will finish the 
job: It provides a way to share the momentum trans- 
fer in such a way that a slow proton recoils while two 
fast nucleons are emitted in the forward direction, with a 
small enough relative momentum to recombine into the 
deuteron. This study remains to be done. 

In the same experiment, the np continuum has been 
recorded. Two body short range correlations are the pri- 
mary source of high momentum components in the nu- 
clear wave function. They are strongly coupled to high 
energy states in the continuum, where they induce a peak 
(dot-dashed line in Fig. [HJ characteristic of the disinte- 
gration of a NN pair at rest in '^He 1]. The width of 
the peak reflects the Fermi motion of the pair. Again, 
FSI between the two nucleons of the pair (dotted line) 
dominate the cross section. The subsequent scattering 
of one of these nucleons with the spectator third nucleon 
(dashed line) shifts the peak toward the experiment 28j , 
but is not dominant. MEG and A formation (full line) 
brings down the cross section in good agreement with the 
experiment. In the continuum, one measures the transi- 
tion between a correlated pair in the "^He ground state 
and a correlated pair in the continuum. It turns out that 
pn pairs (in T = and T — \ isospin states) as well as 
pp pairs contribute by roughly the same amount. 
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FIG. 6: The recoil energy distribution in the ^He(e,e'p)np 

reaction at p™= 622 MeV/c, X=l and = 1.55 GeV^ 



The full implementation of the actual SAID amplitudes 
poses the problem of their extrapolation in the unphysi- 
cal region (off-shell NN scattering) which is under study. 




Triple coincidence studies |29|, |30j of the reaction 
■^IIe(e,e'NN)N have been completed in the full phase 
space with the large acceptance spectrometer CLAS at 
JLab. The model reproduces also the cross section for 
various cuts in the phase space. I refer to my talk at the 
Lisboa conference !31|] for a comparison with preliminary 
data. 

Finally, the model gives a good account of the cross 
section of the reaction '*He(e,e'p)T, which has been re- 
cently determined at JLab: see ref. [33j for a comparison 
with preliminary data. 

It is remarkable that the cross sections of so many 
channels are reproduced with the simple choice (eq. (JJ 
of the central part of the nucleon-nucleon scattering am- 
plitude. While its value at the very forward angle is 
fixed by unitarity, the slope parameter (3 has been de- 
termined by fitting the angular distribution of the unpo- 
larized nucleon-nucleon scattering cross section. It turns 
out that this form reproduces fairly well the modulus 
of the central part of the NN amplitude extracted form 
the SAID data base js^l up to four momentum transfer 
-t = 0.4 GeV^, i.e. pm ~ V-i = 0.63 GeV/c in the 
quasi-elastic kinematics. However, the SAID ratio e be- 
tween the real part and the imaginary part of the scatter- 
ing amplitude varies from about zero at pm = to about 
one around Pm = 200 MeV/c and back to zero in the 
range 300 < Pm < 700 MeV/c. Thus, the parameteriza- 
tion 13 is very good in the recoil momentum range where 
FSI dominate. Above, A formation and MEG take over, 
and the details of the NN amplitude are less important. 



FIG. 7: The ratio of the spin transfer coefficients (top) and the 
induced proton polarization (bottom) in *He(e,e'p)T. Dashed 
lines: PW. Dash-dotted lines: FSI. Dotted lines: 2 body MEG 
and A included. Full lines: 3 body mechanisms included. 



While the spin-orbit and spin-spin terms are taken into 
account in the phase shift expansion (low energy) they 
have not yet been implemented in the absorptive ampli- 
tude (high energy). The spin-orbit term has been found 
to give small contribution to unpolarized observables in 
the (e,e'p) channels [s^l- The present version of the 
model predicts deviations from PW of the ratio of spi n 
transfer coefficients (as defined in Appendix A of llSl) 
of a few % in the direction of experimental values |35j 
recently recorded at JLab (Fig. [TJ . The good agreement 
with the induced polarization Py gives confidence on the 
treatment of various interaction effects. However, before 
drawing any definite conclusion, one has to wait until 
the full implementation of spin dependent terms as well 
as the averaging over the experimental acceptance. 

To summarize, a fair agreement with the recent JLab 
data has been reached around A" = 1, up to recoil mo- 
mentum of the order of 1 GeV/c, provided that the NN 
scattering amplitude relevant to the same energy range 
as well as realistic few body wave functions are used. The 
perpendicular kinematics offers a robust starting point to 
study the evolution with of the re-interaction of nu- 
cleons, but also of hadrons, in view of determining their 
structure at short distances It is not the right 

place to determine the high momentum components of 
the nuclear wave function. One has to go away from 
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the quasi-elastic kinematics: as demonstrated in Fig. ^ 
this occurs in paraUel or anti-parallel kinematics, where 
on-shell nucleon rescattering is suppressed. 

I acknowledge the warm hospitality of Jefferson Lab- 
oratory, where this work was completed, as well as the 
numerous discussions which I had with D. Higinbotham, 
K. Egiyan and E. Voutier over the past few years and 
which have shaped this project. 
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